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Работа посвящена описанию метода адаптивного повышения контраста цветных изображений 
на основе применения сингулярного разложения к результатам нечеткой кластеризации в композиции 
с исходными данными. Предложенный метод позволяет обеспечить сохранение исходной 
цветопередачи снимков. Формирование конечного результата производится на основе автоматизации 
анализа и отбора наиболее значимых составляющих матрицы левых сингулярных векторов. 
Представлены алгоритм и результаты экспериментальных исследований на примере числовой модели 
и реального дерматоскопического медицинского изображения. 
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Введение 

Необходимость повышения качества цветных изображений, в частности, их 
контраста, возникает во многих практических задачах, например, в медицине, 
спутниковом мониторинге земной поверхности, металлографии и пр. Наиболее 
сложной проблемой являются случаи, когда исходное изображение содержит 
участки с различным диапазоном яркостей (контрастные и низкоконтрастные), а 
области интереса (аномалии) несущественно отличаются от фона или шума. Кроме 
того, априори неизвестно наличие и месторасположение аномалий, которые могут 
находится в различных областях цветовой палитры, причем цветовая составляющая 
может быть определяющей при диагностике/анализе. В таком случае, сохранение 
исходного цветового соотношения является принципиально важным. Множество 
существующих методов повышения контраста направлены на решение частных 
задач [1], их применение обычно приводит к искажению исходного распределения 
яркостей и, как следствие, цветопередачи, поэтому разработка новых (модификация 
существующих) методов повышения контраста с сохранением цветопередачи 
остается актуальной задачей. 

В настоящее время одним из эффективных подходов для повышения качества 
изображений является их обработка в пространстве нечетких функций 
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принадлежности [2]. В работе [3] процедура перехода в нечеткое пространство 
(фаззификация) осуществляется на основе метода нечеткой кластеризации. В таком 
случае, получаемый результат, с одной стороны, существенным образом зависит от 
применяемого метода кластеризации, а с другой, многомерность используемой 
нечеткой функции принадлежности, отсутствие универсального метода определения 
значимости отдельных кластеров и неоднозначность выполнения процедуры 
дефаззификации, а, соответственно, и визуализации [4] являются ключевыми 
проблемами, решение которых определяет информативность конечного результата. 

Постановка проблемы 

В работе [5] был предложен алгоритм дефаззификации с использованием 
ортогонализации методом сингулярного разложения. Процедура ортогонализации 
позволяет разделить информацию, что, с принципиальной точки зрения, открывает 
дополнительные возможности повышения чувствительности при обработке 
изображений, устранять избыточность и повышать информативность исходных 
данных за счет приведения вещественной матрицы к каноническому виду, 
отражающему ее геометрическую структуру [6]. 

Суть алгоритма заключалась в том, что после осуществления сингулярного 
разложения для полученной нечеткой функции принадлежности производилась 
автоматическая оценка информативности столбцов матрицы левых сингулярных 
векторов на основе матрицы правых сингулярных векторов. В работе [5] за счет 
слияния результатов локальных (для окон преобразования) и глобальных 
сингулярных разложений было достигнуто повышение чувствительности 
сегментации полутоновых изображений. Однако, непосредственное применение 
этого подхода к визуализации результатов кластеризации цветных снимков 
приводило к нарушению цветового баланса. 

В данной работе предложен адаптивный метод повышения контраста цветных 
изображений на основе слияния глобального и оконного сингулярных разложений 
цветного изображения, полученных в результате нечеткой кластеризации, с 
автоматическим отбором наиболее значимых столбцов матрицы левых сингулярных 
векторов и композицией с исходным изображением, который позволяет сохранить 
цветовой баланс. 

Изложение основного материала 

Предложенный алгоритм повышения контраста цветных изображений состоит 
из следующих шагов: 

1. Выполняется нечеткая кластеризация исходного цветного изображения / 


размерности 4ухахх& (4у и Ах - размерность снимка по вертикали и 
горизонтали, соответственно, (Ё - количество цветовых/спектральных 
компонент) в результате чего формируется нечеткая функция принадлежности И 
размерности ихс (п=ау-Ах, с - количество нечетких кластеров). 


Производимое на этом шагу расширение пространства входных признаков 
позволяет учитывать дополнительную информацию при сегментации. 


2. Формируется матрица /[ с размерностью пх(& + с), содержащая матрицы / и 0. 


1 ы 
3. Для матрицы Г выполняются следующие деиствия: 
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а) 


6) 


в) 


г) 
д) 


е) 


3) 


и) 


к) 
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Осуществляется сингулярное разложение [6], в результате чего формируются 


матрицы /1 т У; (левые и правые сингулярные вектора, соответственно), $5; 
(собственные числа): 

п=П. 5,9). (р 
При этом размерность у матриц /[ Ти / . совпадает, а для матриц 5; и \; - 
имеет значение (а +с)х (4 +с) (используется экономное представление 
сингулярного разложения). Преобразование (1) обеспечивает переход к 
новому пространству информативных признаков. 
Все столбцы [ 1 масштабируются на отрезок [01] и интерпретируются как 
матрица размерности фухахх (@ + с) : 
На основе матрицы У, вычисляется вектор коэффициентов С по следующей 
формуле [4]: 

+ + 


С; =| | У. (и, 1+1 У. (и, п 2+ с]. (2) 
УЕ 1= 


Вектор коэффициентов С упорядочивается по убыванию, а его элементы 
нормируется так, чтобы их сумма была равна 1. 
Формируется вектор АС, содержащий разности для каждой пары соседних 
элементов вектора С. 
Вычисляется значение 4С„ по формуле: 
-+с-1 
ас. =|| У4С,/(+с-1|+(@Сы +АСнах )/2 |/2, (3) 
Я 

где АСпа и АСпих -— минимальный и максимальный элементы вектора 4С, 
соответственно. 
На основе значения 4С, при просмотре элементов вектора 4С с конца 
выбирается такой индекс ах ‚ для которого выполняется условие: 


ас.  >ас.. (4) 


Элементы вектора АС с индексами от 1 до пих +1 нормируется так, чтобы 


тах 


их сумма была равна 1. 


Формируется матрица 1 2 размерностью 4ухахх4 по следующим 
формулам: 


пах 
(2 => (п ) (ас) у = ау х = [Ьах в [2,ай, (5) 
о = № р Е [, ау х = [, 4] 7 = 2,41. (6) 


Для учета влияния последнего столбца матрицы левых сингулярных 
векторов, который, как показали эксперименты, может содержать 
информацию об аномальных объектах (участках) на изображении, 


выполняется следующее преобразование матрицы / ^. 
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2 2 1 
( р = ( а -(-Кеа)+ [п и Копа» УЕ [4х4], (1) 
Кена =(@С> +4С; „. +АСцье 3. (8) 
Матрица 1 ь интерпретируется как цветное изображение, каждый цветовой 
канал которого масштабируется на [0,1]. 


4. Матрица 1 т разбивается на неперекрывающиеся окна заданного размера, и для 
каждого окна выполняются следующие действия: 
а) применяется сингулярное разложение, в результате чего вычисляются 


1 
матрицы \Г5 у’иа ? (у, а ь (5, а 
1 1 Т 
( лы -(' 5 Пе -(5 5 д - (у, д : (9) 
При этом размерность у ( } И п „а СовВпадает и соответствует 


Ата хаха х(@+с), где 4уиа и Ча - количество пикселей у 
текущего окна по вертикали и горизонтали, соответственно. Размерность 
матриц ($, а ай а 4 Имеет значение (а +с)х (4 +с). 


6) Для матрицы ( . 4 повторяются все действия из подпунктов 6) — к) шага 3. 
„ 13 
В результате слияния матриц оконных преобразований [ Я формируется 


матрица /[ 3 размерностью 4у.а4хх а. Эта матрица также интерпретируется как 
цветное изображение, каждый цветовой канал которого масштабируется на [0,1]. 


При выполнении шага 4 вычисляется матрица АС„и4, которая является 
усреднением всех полученных матриц АС„иа. 
5. Формируется матрица 1 . путем слияния матриц, полученных в результате 


применения «оконного» и «общего» сингулярных разложений, что выполняется 
для подавления «оконного» эффекта: 


ПЕ к: К+ р ры -1- К), уе, 4 хехе ай, (10) 
Рух = О Че + 3: Кута, УЕ [ау] х Е 1.4]: Е [,&&], (11) 


причем коэффициенты К; 


т? И Купа (для каждого окна преобразования) 


вычисляются По одной и той же методике, но с использованием разных векторов 


АС (для вычисления К;„, используется вектор АС, полученный при 


формировании матрицы 1 т: а для К,„иа - вектора АС текущего оконного 
преобразования) согласно следующих формул: 


ЕК. 
Кипа /1те = (К. ук.) , (12) 
Кс 
Ки = пи (К, к, к. нас) } (13) 
АС1 + аС; 
К. = К; .( 1 т (14) 


2 


Илах 
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пах 


К; = >. аС,, (15) 

1=1 

а коэффициент К рассчитывается так: 
К=0.25+(Ко) и Кс, (16) 


причем (Кеть вычисляется по формуле (14) на основании вектора 4С для 


матрицы 1 2. а К. — среднее значение коэффициентов Кс, полученных по 
формуле (14) для каждого окна преобразования. 


6. Цветовые каналы выходного изображения /[ 2 формируются путем слияния 


цветовых каналов матриц /[ Чи Л шо следующим формулам: 


т +11 (-Кв), УЕ [4 хе], (17) 
=К 14 
Кв =(к! | т (18) 
РИ > ЕЕ 
к = ас + > ас, к) 12 (19) 
1=4-+1 


для красного цветового канала. При этом й„.х - усредненное (с округлением) 
значения индекса [„х для оконных преобразований. Для зеленого и синего 


цветовых каналов формулы преобразования следующие: 
4 4 
т = м +Гух: К; /2+1,х;а-К, уе ау] хе [вах [2,4!|, (20) 


к; = (ас; + Кв ий Е [2,4 . р 

Экспериментальные результаты 

Предложенный алгоритм был применен при обработке различных цветных 
изображений, примером которых являются снимки, представленные на рис. 1, 2. 

На рис. 1 а приведено цветное модельное изображение, содержащее на 
неравномерном фоне 7 визуально неразличимых аномальных объектов 
(геометрические фигуры). На рис. 1 в -— 1 г приведены цветовые каналы этого 
изображения, анализ которых показывает, что аномалии отчасти видимы только на 
синем цветовом канале (рис 1 г), но неразличимы на цветном изображении. 

На рис. 2 а приведено медицинское изображение, представляющее собой 
результаты дерматоскопии, проведенной с целью диагностирования меланомы по 
визуальным признакам, среди которых основными являются наличие бело-синих 
структур и пятен неправильной формы. На рис. 2 6 - 2 г представлены цветовые 
каналы этого изображения, а на рис. За-— Зв- их гистограммы. 

При проведении экспериментов был использован алгоритм нечеткой кластеризации 
$НСМ [6, 7]. При его применении динамическое сжатие и обработка нечеткой функции 
принадлежности не использовались. Перед применением карты Кохонена число нечетких 
кластеров увеличивалось в 4 раза, а также осуществлялось упорядочение центроидов на 
основе матрицы Евклидовых расстояний. Выбор исходного числа нейронов после 
обучения карты Кохонена осуществлялся на основе максимального соответствия. 
Количество нечетких кластеров с =6 в обоих случаях. Размер окна при применении 
предложенного в данной статье метода был выбран 4х4 пикселя. 
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а б 
в г 


Рис. 1. Модельное цветное изображение (200х200): а — исходный снимок; 
цветовые каналы: 6 — красный, в — зеленый, г — синий. 


В Г 
Рис. 2. Цветное медицинское изображение (323х215): а — исходный снимок; 


цветовые каналы: а — красный, б — зеленый, в — синий 
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Рис. 3. Гистограммы цветовых каналов изображения рис. 2а: 
а - красного; б — зеленого; в — синего 
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На рис. 4 а представлено цветное изображение, полученное после повышения 
контраста предложенным методом у снимка на рис. 1 а. При сохранении цветового 
баланса по сравнению с исходным снимком было также достигнуто четкое 
визуальное выделение 3 скрытых объектов в левой части снимка и визуальную 
различимость 4 скрытых объектов в правой части. При этом визуализация 
результатов нечеткой кластеризации на основе сравнения с исходными данными [8] 
(рис. 4 д — цветное изображение) не позволяет выделить «аномальные» объекты, а 
также приводит к нарушению цветового баланса. 


г д 
Рис. 4. Повышение контраста изображения рис. Та: а — результирующее цветное 
изображение; его цветовые каналы: 6 — красный, в — зеленый, г — синий; 
д - визуализация на основе сравнения с исходными данными (цветное изображение) 


Применение предложенного метода к изображению на рис. 2 а (рис. 5 а) 
обеспечивает повышение контраста в области интереса  (обведена 
прямоугольником), что позволяет диагностировать наличие меланомы на основании 
присутствия бело-синих структур. При этом сохраняется цветовой баланс, что 
можно увидеть при сравнении цветовых каналов полученного изображения (рис. 5 6 
—5 ги их гистограмм, приведенных на рис. ба - 6 в, с цветовыми каналами (рис. 2 
6 — 2 г) и их гистограммами (рис. 3 а - 3 в) у исходного снимка. Визуализация 
результатов нечеткой кластеризации на основе сравнения с исходными данными 
(рис. 5 г), хоть и позволяет обнаружить наличие бело-синих структур, но приводит к 
нарушению цветового баланса, что затрудняет визуальное определение области 
распространения меланомы. 

Следует отметить, что при применении предложенного метода с целью 
уменьшения различимости оконного эффекта следует выбирать небольшой размер 
окон. Рекомендуемый диапазон: от 3х3 пикселя до 6х 6 пикселей. 
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Рис. 5. Повышение контраста изображения рис. 2а: а — результирующее цветное 
изображение; его цветовые каналы 6 — красный, в — зеленый, г — синий; 
д - визуализация на основе сравнения с исходными данными (цветное изображение) 
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Рис. 6. Гистограммы цветовых каналов изображения рис. 5а: 


а — красного; б — зеленого; в — синего 
Выводы 
Предложенный в данной работе адаптивный метод позволяет при визуализации 
результатов нечеткой кластеризации цветных изображений как добиться повышения 
контраста (в некоторых случаях также достигается повышение яркости), так и 
сохранить цветовой баланс по сравнению с исходным снимком. Достоинствами 
предложенного метода является адаптивность и отсутствие управляющих 
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параметров. К недостаткам следует отнести недостаточный уровень подавления 
оконного эффекта, что проявляется при визуальном анализе относительно 
небольших изображений. 
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